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We present a set of modified quantum rate equations, with the help of the nonequilibrium Green's 
function and slave-particle techniques along with the correct quantization, for description of the 
quantum transport through an interacting mesoscopic region connected with two leads, in the se- 
quential tunneling regime. The assumption that only leading order of \V\^ {V is the tunneling 
coupling between the interacting central region and the leads) has been taken into account in de- 
riving these equations implies that the quantum rate equations are only valid in the case of weak 
coupling between the central region and the leads. For demonstrations, we consider two special 
cases in the central region, a single interacting quantum dot (SQD) with weak spin-fiip scattering 
and a weakly coupled double quantum dots (CQD), as examples. In the limit of zero temperature 
and large bias voltage, the resulting equations are identical to the previous results derived from 
the many-body Schrodinger equation. The numerical simulations reveal: 1) the dependence of the 
spin-flip scattering on the temperature and bias voltage in the SQD; and 2) the possible negative dif- 
ferential conductance and negative tunnel magnetoresistance in the CQD, depending on the hopping 
between the two quantum dots. 

PACS numbers: 73.21. La, 73.23.-b, 73.23.Hk, 85. 35. Be 



I. INTRODUCTION 



For many years, much experimental and theoretical 
work has been devoted to exploring the transport prop- 
erties of artificially nanofabricated structures containing 
a discrete number of quantum states and a small num- 
ber of electrons. The tunneling current through these 
mesoscopic devices, isolated from two macroscopic leads 
by potential barriers, manifested many novel effects due 
to this confinement. For example, in a semiconductor 
quantum dot (QD) one observed the Coulomb blockade 
oscillations due to the charging energy^ and even the 
Kondo effect due to the strong on-site Coulomb inter- 
action in the tunneling transport. ^'"^ Recently, interest 
in quantum computation and quantum information pro- 
cessing has attracted increasing attention to the problem 
of measurement of tunneling currents via a mesoscopic 
system that can be modeled by a two-level Hamiltonian, 
for example, charges in coupled QDsi^i^ and spins in 
a QD under magnetic fieldsi^ Measurement of the tun- 
neling current in such systems provides information not 
only about the Rabi oscillations^ between the two lev- 
els, but also about the spin precession in quantum spin 
oscillationS(2iiS both of which are crucial improvement in 
the control of the superposition of the quantum states 
and thus quantum information processing. In addition, 
similar physical picture has been utilized with success 



to analyze transport through molecular nanojuction,li 
for example, a system of benzenei^ and DNA molecular 
chain li^ 

In order to describe this kind of quantum oscillations in 
quantum transport through a QD, the master equations 
or "quantum" version of the rate equations have been 
first proposed by Nazarov,''-^ and later derived microscop- 
ically from the Schrodinger equation directly^^*^ and 
from the von Neumann equation and superoperatorsii, 
respectively. In these, the central point is that the equa- 
tions of motion (EOM) of the diagonal density matrix 
elements allow for an additional term of the nondiag- 
onal density matrix elements, which indeed stand for 
the coherent superposition of different quantum states 
and are referred to as the coherent transfer term, along 
with the time evolutions commanded by their own EOM. 
We have to solve these equations selfconsistently to de- 
termine the nonequilibrium probability densities. As a 
result, the tunneling current unavoidably contains the 
contributions of the nondiagonal density matrix elements 
and naturally provides the information of the quantum 
Rabi oscillation, although the explicit expression of cur- 
rent formulation only involves the diagonal density ma- 
trix elements. The modified quantum rate equations have 
proved successful in describing this quantum oscillation 
in coherently coupled quantum dots (CQD)fi4 the quan- 
tum measurement by using quantum point contact near 
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QQp^i8ji9j20 gj^^ even time-dependent quantum tunnel- 
ing through the CQD.^^ On the other hand, the Coulomb 
interaction inside the small confined region plays cru- 
cial role in, as above mentioned, determining the quan- 
tum transport properties of the devices, of course in con- 
trolling the quantum oscillations of two-level systems. 
In fact, the so-called noninvasive quantum measurement 
process is also based on the Coulomb coupling between 
the detector and the measured systemiSi^ To our knowl- 
edge, however, a systematic investigation of the quantum 
rate equations at arbitrary temperature associated with 
the Coulomb interactions has been lacking. 

About ten years ago, a "classical" rate equation 
was derived for sequential tunneling through a double- 
barrier system from the "quantum" kinetic equation, the 
nonequilibrium Green's function (NGF), which is be- 
lieved to be a more powerful tool for studying nonequi- 
librium phenomena. ^^'^ ^ Our aim in this paper is to sys- 
tematically explore the quantum rate equations for the 
interacting systems in the sequential tunneling regime 
from the NGF approach. The unique source of dif- 
ficulty is how to deal with the Coulomb interaction 
term in the derivation. This problem is the same 
as that we have in studies of the strongly correlated 
fermionic systems, for example, the recent investigations 
on the Kondo enhanced conductance of a QD at low 
temperatureiSiSSiSi' Many theoretical methods have 
been developed to solve the strong correlation effects. 
Among these methods, the slave-particle technique is of 
particular elegancei^SiSLS^i'^ The great advantage of this 
approach is that the correlated Hamiltonian for the sys- 
tem under studied is transformed to an equivalent one 
without Coulomb correlations, while introducing several 
auxiliary particles. Thus previously well-developed for- 
mulations for noninteracting systems can be applied to 
investigate the interacting systems in the framework of 
the new representation. Along this line, a further tech- 
nique advance is made in the present work. Here we 
extend the approach of the slave-particle representation 
to the weakly coupled quantum system of interest here 
and give the consistent Hamiltonian formulation in terms 
of the slave particles. The equations of motion of the 
density matrix elements are then studied in the frame- 
work of NGF and within the slave-particle scheme. Our 
derivation contains three approximations. The first is to 
assume that the central region has very "weak" coupling 
with the external environments (the leads) V. Secondly, 
we assume the couplings between subsystems are also 
weak to keep them being individual, for example, the 
weak spin-fiip scattering in SQD and weak interdot hop- 
ping in CQD. As a result, we can give the definitions 
of the spectrum expressions of the NGFs of the central 
region in terms of the nonequilibrium probability densi- 
ties and keep only the leading order term in jyj^ in the 
expansions of the equations of motion. The final one is 
to apply the wide band limit in the two leads, namely 
that the coupling strengths between the central region 
and the leads are independent of the energy and can be 



considered as constant. 

The rest of the paper is organized as follows. In the 
next section, we give the derivations in detail for a sin- 
gle interacting QD (SQD) taking the weak spin-flip scat- 
tering into account, and establish the temperature and 
bias voltage dependent quantum rate equations for arbi- 
trary Coulomb interaction. In section HI we derive the 
quantum rate equations for the weakly coupled QDs. In 
both of the two section, after the analytical results are 
discussed for the no doubly occupied level and the deep 
level situations, we perform numerical simulations on the 
occupation numbers and the tunneling current in the gen- 
eral case as functions of the bare level in the QD and the 
bias voltage between the source and the drain. Finally, 
all the results are summarized in Section V. 



II. SINGLE QUANTUM DOT 

We begin with our derivation of the quantum rate 
equations for a SQD with a weak spin-flip scattering in 
this section. In the case of no spin-flip terms, the rate 
equations are of "classical" variety, which have been ad- 
equately described by other methods. Therefore our re- 
sults are not new in this case but are established from a 
different scheme. The purpose of this section is also to 
provide an examination to prove this approach in com- 
parison with previous results in no spin-flip case. 

A. Model Hamiltonian and slave-particle 
representation 

We use the standard model Hamiltonian to describe 
the resonant tunneling through a SQD, as depicted in 
Fig. 1(a), with a single bare- level ed and a weak intradot 
spin-flip scattering connected to two non-interacting 
leads: 

+ Vui^na^ + ^ ^^oc\^,„c^^ + H.c), 

(1) 

where cj^j,^ {crika) and c^^ icda) are the creation (annihi- 
lation) operators for electrons with momentum fc, spin-cr 
and energy trtka hi the lead r\ (= L, R) and for a spin-cr 
electron on the QD, respectively. The third term de- 
scribes the Coulomb interaction among electrons on the 
QD. ~ (i.a'^da occupation operator in the SQD. 

The fourth term represents the tunneling coupling be- 
tween the QD and the reservoirs. We assume that the 
coupling strength V^a is spin-dependent, being able to 
describe the ferromagnetic leads. Each of the two leads 
is separately in thermal equilibrium with the chemical 
potential /x^, which is assumed to be zero in equilibrium 
condition and chosen as the energy reference throughout 
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the paper. In the nonequihbrium case, the chemical po- 
tentials of the leads differ by the applied bias. In this 
paper, we assume the tunneling coupling is weak enough 
to guarantee no Kondo effect in our model and the QD 
is in the Coulomb blockade regime. Generally, we have 
_R| = i?J = i? being a constant. 

According to the finite-?/ slave-boson approach intro- 
duced by Zou and Anderson^ the electron operator Cda 
can be written in four possible single electron states, 
namely: the empty state |0) with zero energy eg — 0; 
the singly occupied (with spin up or down) electronic 
state \a) with energy = ^d, and the doubly occupied 
state I ti) with energy Sd — 2ed + U, as 



Cd. = \0){<j\+a\a)(n\, [c 



(2) 



Because these four states expand the entire Hilbert space, 
the completeness relation must be satisfied 



|o)(oi + 1 nxn I + Ei^)H = i- 



(3) 



These Dirac brackets were then treated as operators: 
e'' = |0), = I ID as slave-boson operators and = |(t) 
as pseudo-fermion operator. In terms of these auxiliary 
operators, Eqs.(|2l and Q become 



(4) 
(5) 



The explicit (anti)communicators of these auxiliary par- 
ticles can be easily established from the definitions of the 
Dirac bracketsi^ 



ee' = 1, dd^ = 1, f^fl, = Saa', 



J - 
ed^ 



efl = Le^ - Ld^ - det ^ rf/t = 0. (6) 



Therefore, along with these correct quantization, the 
Hamiltonian ^ can be replaced by the following form 
in the auxiliary particle representation: 

+(2ed V)Sd + i?T/|/i + Rifj^ 

+ E + ojU) + H.c], (7) 

which was proved to be equivalent to the original one 
by Zou and Anderson in the case of no spin-flip termiSS 
Furthermore, as far as the four possible single electric 
states are considered as the basis, the statistical expec- 
tations of the diagonal elements of the density matrix. 
Pa [i — {0,cr, d}), give the occupation probabilities of 
the resonant level in the QD being empty, or singly occu- 
pied by spin-cr electron, or doubly occupied, respectively. 
The nondiagonal term describes the coherent super- 
position state between the spin-up and -down states in 
the QD. In the slave particle notation, the correspond- 
ing relations between the density matrix elements and 



these auxiliary operators are obvious as poo = 
ete, p^^ = \a){a\ = /J/,, pdd = \ \ = Sd ' and 

the nondiagonal term p„s = |ct)((t| = flfa- According to 
Eq. (O, the constraint is subject to the diagonal elements 
of the density matrix poo + P<y<y + Pdd = '^- 



B. Derivation of the quantum rate equations 

In this subsection, we derive the rate equations for se- 
quential tunneling starting from the combined fermion- 
boson Hamiltonian |7J by using the Keldysh's NGF. 

In order to describe the nonequihbrium state of 
electrons, we define the retarded (advanced) and 
lesser (greater) Green's functions (GFs) for the QD 

((c,.(t)|cL{i')))'^(")'<(>) as follows: 



G:ii{t,t') EE ±^e{±tTt'){{c,^{t),cl^,it')}), G<At.t') EE 

X4'(<')c,.W) and G>„,{t,t') EE -^{c,^{t)cl,{t')). 
Considering Eq. I^J, these GFs in the QD can 
be divided into two parts Ga^' — Geaa' + Gdaa' 
with Ge..' = {{eKt)f,{t)\fl{t')e{t'))) and Gdaa' - 
aa' {{fl{t)d{t)\d'^ {t')f^,{t'))). Under the weak coupling 
assumption, the central region can be regarded as an con- 
siderably isolated system and its density matrix operator 
is supposed to be p = PijPij = {0, o, d}), mean- 
while the reservoirs are taken as "environment" located 
in local thermal equilibrium. Therefore, we can readily 
define the decoupled diagonal GFs of the QD for weak 
spin-flip transitions in terms of spectrum expression, in 
the Fourier space, as 



Gllai^) = 



Poo 



G<a'a{^) 

G^U^) 

G^Jai^) 

G,>°,(^) 



Lo — ed + *0+ ' 

Pdd + Paa 



uo-{ed + U)+ iO+ ' 
= 2mp„„5{uj - e<j), 
= 2iTipdd5[io - {ed + U)], 
= -2mpoo5{uj - ed), 

= -2TTipsg5[uJ - {ed + U)]. 



(8) 



If no bias voltage is added between the two leads, the 
central region is in a thermal equilibrium state, and the 
distribution probabilities are well-known as poo = 1/^, 
p,, = e-'^/^/Z, and pdd = 6-^^'^+"^^ /Z with Z = 
1 + 2e^^'^^^ -I- e~(2'^<i+^)/'^. As far as the spin-flip transi- 
tion is considered, the decoupled nondiagonal correlation 
GFs are crucial in the following derivation. Their Fourier 
expressions are easily obtained from the definitions as: 



G'^aA^) = 27ripCT5(5(w - e<j), 
G>o,(c.) = -2^ip,-,8{^-{ed^U)l 
G<Uu;) ^ 0, G>,%{u;)^0. 



(9) 



For the case of nonequihbrium interested here, these out- 
of-equilibrium probabilities are determined by the cou- 
pling to environments with different chemical potentials, 
and usually they obey a set of equations of time evolu- 
tion, the rate equations. 
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Here, we start from the equations of motion of the op- 
erators pij with the Hamiltonian Q and modified quan- 
tization Eq. JH)): 

Poo = i[H,e^e] = -i^^ (^rj^r^^^^eV^ - y^^flec^ka)^ 

(10a) 

7],k 

+iRsflf, - iRaPj^, (10b) 

71, k 

+iRAfU.~flh), (lOc) 

I 



Paa 



Pas 



in terms of the QD's GFs in the Fourier spacef^^ Here 
r,,cr = 27r^^ |T^crP^(w — €r]kcr) denotes the strength of 
coupling between the QD level and the lead "q. In wide 
band limit, it is independent of energy and is supposed 
to be constant. Under the weak coupling assumption, it 
is adequate to keep only the leading order of 1^1^ in eval- 
uation of these occupation densities. So we can replace 
these QD's GFs with their decoupled formulas Eqs. ^ 
and Q. Finally, the resulting quantum rate equations 
become: 

Poo = ^{T^Paa - r+poo), (12a) 



Paa = r+poo + Tff Pdd ^ a + ^t)P<Ja ^ '^.^{RsPcrs), 

(12b) 



Pdd = i[H, d}d] = i ^ ^{Vr^aclkafU - V^ad'^ fsCrika)- 
rj.k,(T 

(lOd) 

Their statistical expectations involve the 
time-diagonal parts of the correlation func- 
tions: G<^,,,.., (t, i') ^ ^(^-^ (i')e^ W/. W), 
Gt>,,k.it^t') ^ ^{cl^{tVlmt)), G<^,,,(t,t') ^ 
^{n{t'Ht')c„,,,{t)), and G<,^^,^,{t,t') 

^{d^{t')fa'it')Cnkait))- With the help of the Lan- 
greth analytic continuation rules(2i we obtain the 
following expressions in the wide band limit (The detail 
derivation will be given in the Appendix): 



r 



p^g = ii?^(PCTCT-Pffff)--(r++rJ+r„ +r^)p^5, (i2c) 



Pdd = rfpTT + Pu - (^T + -'"i ^P'i'i' (^^^) 

together with the normalization relation poo + Pdd + 
Paa — 1 from Eq. with the definitions F^ — 
T.r^'^v'^frfi'^d) and f± = E,,r»?<T/^(ed + U), where 
f+{u;) = {1 -I- e(""''")/^}"i is the Fermi distribution 
function of the rj lead and /^(w) = 1 — f~^{oj). Here, 
r+ (F^) describes the tunneling rate of electrons with 
spin-CT into (out from) the QD without the occupation of 
the (T state. Similarly, F+ (F") describes the tunneling 



J d^^{F,,/,(^)G>,(c.) +F,^[1 - /,(c.)]G<,(c.)}, (11a) 



P"" = 2^ 



^ J rf^^{F,./,(^)G> ,(c.) + F,.[l - f,iLo)]G<,,{Lo) - F,5/,(c.)G>,,(c.) - F,^[l - Uio)]G<^jLu)} 

+iRaPaa - iRaPaa, (Hb) 

= f / '^'^E(V + r,„^){/,(c.)G> ^(c.) + [1 - f,iuj)]G<,^iuj) - f,iuj)G>M - [1 - ./,(c^)]G<,,(c^)} 

+iRaiPaa - Paa), (Hc) 



Pdd^^J diuJ2{^v'rM^)G^sA^) + r,,[i - /,(^)]G<^J, (lid) 



r\.a 
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rate of electrons with spin-cr in to (out from) the QD, 
when the QD is aheady occupied by an electron with 
spin-CT, revealing the modification of the corresponding 
rates due to the Coulomb repulsion. 

These rate equations Eqs. (|12a|l . Ijl2b|l and Ijl2d|l coin- 
cide with the previous classical rate equations in the se- 
quential picture for the resonant tunneling if the intradot 
spin-flip transition is quenchediSi^ Obviously, if the left 
lead has the same chemical potential as the right lead, 
the stationary solutions of Eqs. p2a|l . (|12b(l and l|12d(l re- 
duce exactly to the above-mentioned thermal equilibrium 
results in the case of i? = 0. In this situation, they have 
clear classical interpretations. For example, the rate of 
change of the number of the spin-cr electrons in the 
SQD, described by Eq. I|12b(l . is contributed from the 
following four single-particle tunneling processes: 1) tun- 
neling into the QD with spin-tr electrons r+ from both 
left and right leads if the QD is initially in the empty 
state Poo; 2) tunneling out from the QD with spin-a elec- 
trons into both two leads if the QD is initially in the 
doubly occupied state pdd', 3) tunneling into the QD with 
spin-(T electrons from both two leads; and 4) tunnel- 
ing out from the QD with spin-cr electrons T~ into both 
two leads, when the QD is initially just in this state po-o-- 
Tunneling events 1) and 2) increase the probability of 
the spin-cr state, but events 3) and 4) decrease this prob- 
ability. These contributions constitute the classical rate 
equation form. Other diagonal equations have similar 
interpretations. Notice that the final term in Eq. Ijl2bp 
describes transitions between isolated states through the 
coupling with nondiagonal terms, which has no classi- 
cal counterpart. Therefore, it is responsible for coherent 
effects in the transport. 

The nondiagonal matrix element p^g is ruled by 
Eq. H12c|l . which resembles the optical Bloch equation 
and describes the dynamics of quantum superposition. 
This is a pure quantum effect. As mentioned by Gurvitz 
and Pragerfi^ the couplings with the leads (all possi- 
ble tunneling processes involved) always provide negative 
contribution and cause damping of the quantum super- 
position. 

The particle current 1^ flowing from the lead rj to the 
QD can be evaluated from the rate of time change of the 
electron number operator Nr,{t) = Y^k.a ^\kA^)%kA^) 
the lead r] 



.32 



m = 



h\ dt 



k.a 



k.a 



-v,;MtitMt) + <jd^Mt)]c.{t)}). (13) 



Ultimately, the current can be expressed in terms of the 
CPs in the QD: 



. duj 



V[l - Uu)][G<,,{u;) + G%^]}. (14) 
Under the weak coupling approximation, it becomes 

Iv/(^ = ^(^v<^Pdd + ^^aPcyy - ^tsP'yy - ^^aPoo)- (15) 

(T 

This formulae demonstrates that all possible tunneling 
processes relevant to the lead ry can provide correspond- 
ing contributions to the current of the lead 77 and the 
current is totally determined by the diagonal elements 
of the density matrix of the central region. However, 
the nondiagonal element of the density matrix is coupled 
with diagonal elements in the rate equation Eq. Ijl2b|l . 
and therefore influence the tunneling current indirectly. 

C. Discussions 

The rate equations H12|l may be readily solved under 
stationary condition for arbitrary bias voltages V and 
temperatures T , and consequently the dc current may be 
obtained via Eq. (|15|l . More interestingly, it is useful to 
review the following two special cases in the case of large 
Coulomb repulsion. First, we consider that no doubly 
occupied state is available in the QD, i.e., p^a = 0. In 
this case, we assume that the bare-level of the QD 
is just above the Fermi levels p of the two leads under 
equilibrium condition, meaning ~ 0, F" ~ Si;^'?^- 
Then, in steady state, the quantum rate equations (jl2b|l 
and H12c|l become 

F+poo - r^pcrcr - 2'^{Rgp„s) = 0, (16a) 

iRa{Paa - Paa) - ^(F^T + F^)pcra 0, (16b) 

with Poo + X)(T Pctct = 1- They can be readily solved 

F+F,-+x(r+ + F+) 



PTT 



Pii - 



F+F-+x(F++F+) 



A 



(17a) 
(17b) 



in which A = [(F+ + F-)(F+ + F-) - F+F+] + x(2F+ + 

2F+ + F^ + Fp and x = 4|i?| V(F^ + TJ). The steady 

tunneling current is Ir = eY,a\(^Ra+^Rs+^Ra)P'yy - 
F^^]. For large bias voltage, i.e., F^^ = and F^^ = 0, 
and spin-independent tunneling, the dc current becomes 



2Fi(Ffl + 2x') 



2FL+FH + x'(2 + 4FL/Ffl)' 



(18) 



and x' — 2|i?P/Ffl, which depicts the spin-flip transi- 
tion induced modification for the corresponding formula 
Eq. (3.10) in Ref.[li 

The second case we considered here is the deep level 
in the large-?/ limit: the bare- level is far below the 



6 



Fermi level ^ but + L'^ is slightly above the Fermi level 
in equilibrium condition, implicating that the QD is al- 
ways occupied by electrons. In this situation, we have 
Poo = 0; r~ ~ 0, and r+ ~ J2n^v<^- Different from the 
above case, only singly and doubly occupied states are 
permitted in tunneling processes. Under the stationary 
condition, the quantum rate equations Ijl2b|) and (|12c|) 
reduce to 



(19a) 



iRAPaa - Pas) - ^{T^ + rt)Paa = 0, (19b) 

with pdd + X)cr Pctct = 1- The solutions are 



(F+ + x)(r+ + r-) + x(r+ + F-) ^ 



(F+ + x)(F+ + x)-x2 



Pii 



Pdd = 



(F++x)(F++x)-x2 

(f+ + x)(f+ + f-) + x(r+ + r-) 



Pdd, 

(20a) 

Pdd, 

(20b) 



(F+ + x)r+ + xr[ 
-x)(f+ + f-) + x(f+ 



(F+ + x)r+ + xr^ 
withx 4|i?|V(r++f^ 



(20c) 



J . The steady tunneling current 

is given by Ir = eJ^A^RaPdd - ^rbPctct]- Large bias 
voltage further simplifies the spin-independent current 
as 



iR/e- 



2TlTr + 2x" 



Fi + 2F« + 2x"(H-Ffl/Fi)' 



(21) 



with x" = 2|i?p/fL. This is a modification of 
Eq. (3.11) given by Gurvitz and Prager due to spin-flip 
transitions 

It should be noted that the same two cases are also 
analyzed in Ref. |3^ for an interacting QD with spin- 
flip transitions included. They evaluated the occupation 
numbers from the classical rate equations and utilized a 
spin relaxation time Ts to describe the spin-flip transi- 
tions. Therefore, their results are slightly different from 
ours. For both cases, if we redefine the spin relaxation 
time as I/ts = x, which is now of temperature and bias 
voltage dependence, their results'^^ are the same as ours, 
Eqs. (|17|l and (|20|) of the quantum rate equations. This 
is a clear demonstration of the importance of quantum 
"coherence" . 

It is also worth examining the quantum rate equations 
(|12|l derived here at large bias voltage between the left 
and right leads without the spin-flip transitions. In this 



case, we assume eV ^ T and eV ^> U, so F^^ = Pjj^ = 
and F£^ — F^^ = 0. The quantum rate equations lead: 

Poo = Tt^ipii + Fflxpxi - (Fl| + Flx)/5oo, (22a) 

Paa = ^LaPOO + ^RaPdd — Ra + ^Ls)Pcrcr, (22b) 

Pdd = T^LiPu + TLTPii - (rflT + ^Rl)Pdd- (22c) 

And the current is Ir = eJ2a(^R<^Pdd + ^RaPaa)- At 
zero temperature and spin-independent tunneling, these 
equations indeed resemble the rate equations derived 
from the Schrodinger equation developed by Gurvitz and 
Prager .M 

D. Numerical results 

In this subsection, we perform numerical calculations 
for the spin dependence of the tunneling processes, 
through the SQD connected to two ferromagnetic leads. 
In the following calculations, we consider two magnetic 
configurations, namely, parallel (P) and antiparallel (AP) 
configurations. When the magnetic electrodes are in P 
configuration, we assume that the majority electrons are 
spin-up a =t, and the minority electrons are spin-down 
a =i. We also assume that in the AP configuration the 
magnetization of the right electrode is reversed. 

Therefore, for the identical leads and symmetric bar- 
riers, of interest in the present paper, we further assume 
that the fcrromagnetism of the leads can be accounted for 
by the polarization-dependent couplings F^j = F^j = 
(H-p)Fo, Tli = Tri = (1 -_p)Fo for the P alignment, 
while Tli = ^Rl = (1 + P)^o, Tlx = Tr,^ = (1 - _p)Fo 
for the AP alignment. Fq denotes the tunneling coupling 
between the QD and the leads without internal magne- 
tization, and p (0 < p < 1) stands for the polarization 
strength of the leads. In the wide band limit, Fq is sup- 
posed to be a constant and chosen as unit of energy in 
the following paper. Moreover, we measure energy from 
the Fermi levels of the left and right leads in the equilib- 
rium condition (/i^ = PR. = 0) thereafter. The discrete 
level Ed of the QD can cross the Fermi levels by tuning 
the gate voltage in experiments. Without loss of general- 
ity, we apply the bias voltage V between the source and 
drain symmetrically pL = —pR. = eV/2, and neglect the 
shift of the discrete level caused by this external voltage. 
Because of the symmetry, we will restrict to positive bias 
only, V > 0. 

From Eqs. (|12|) . one can find all the expectation values 
of the density-matrix elements for a given bias V in the 
stationary condition, and thus allow us to calculate the 
tunneling current flowing through the system by employ- 
ing Eq. (|15l) and the nonequilibrium occupation numbers 
Tlx, defined by = p^a + Pdd- 

First we consider no spin-flip scattering processes on 
the QD. Figs. 2(a) and (b) plot the nonequilibrium occu- 
pation numbers as a function of the bare level calculated 
for a small bias y = 1.0 and a large bias V = 10.0, re- 
spectively, in both P (thin lines) and AP (thick lines) 
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configurations. Tiie two spacial bias voltages are chosen 
in order liere to demonstrate the linear response regime 
{V = 1.0) and the strong nonlinear case {V — 10.0), 
respectively. For comparison, we also plot the equilib- 
rium occupation numbers in Fig. 2(b). From these figure, 
we can observe that: 1) The complete Coulomb block- 
ade (charging) effect in equilibrium (the single step in 
Pdd) is partially removed in nonequlibrium, i.e., pdd be- 
comes a multi-step function of the gate voltage; 2) 
has fractional steps in nonequilibrium in contrast to just 
half-interger steps in equilibrium; 3) = in the P 
configuration, whereas ^ in the AP configuration. 
Fig. 2(c) shows the tunneling current calculated for both 
configurations. The current in the P alignment is always 
larger than that in the AP alignment in the whole range 
of the gate voltage. In the linear response regime, the 
current provides the information of the conductance of 
the device: there appear two resonant peaks with equal 
heights when the gate voltage controlled levels and 
ed + U respectively cross the Fermi levels of two leads. 
While in the strong nonequilibrium case, there are three 
steps in the current, which correspond to the steps in 
the occupation numbers, whereas between the steps the 
current is constant. 

Fig. 3 illustrates typical variations of the occupation 
numbers and the current with the bias voltage V for ed — 
1 (the no doubly occupied level) and = — 5 (the deep 
level). In the first (second) case, the first step in 
occurs at the bias, when the Fermi level of the source or 
drain crosses the discrete level ed {ed + U). This means a 
new channel open for tunneling. Consequently, we find a 
step in the current appears at this position. As the bias 
further increases, they all keep constant until the second 
step at a higher voltage corresponding to the case when 
the Fermi level crosses ed + U (ed), which also induces 
a step in pdd- The insets in Fig. 3(c) and (f) depict the 
corresponding tunnel magnetoresistance (TMR), defined 
as 

TMR=^^—^. (23) 
Iap 

The TMR is enhanced by the Coulomb interaction in 
the range between the two biases corresponding to the 
two steps in the current. In these figures, we also dis- 
play the temperature effect in tunneling characteristics. 
It is easily observed that increasing temperature gradu- 
ally smoothes the steep step structure in the occupation 
numbers and the current, and decreases the TMR. 

We now consider the effect of spin-flip scatterings on 
the tunneling. Because the spin-flip processes have no 
influence on the occupation numbers and the current in 
the P configuration, we plot the calculated results for 
the AP configuration with i? = 1 in Fig. 4. It is obvi- 
ous, in comparison with the case of no spin-flip scattering 
R = (thin lines), that the spin-flip transition decreases 
the difference between n-f and n^, increases pdd and the 
current. Moreover, their temperature behaviors are sim- 
ilar with the case of no spin-flip transition. It is worth 



noting that when the bias voltage is lower than 10.0, i.e., 
the value corresponding to the second step in ria- and the 
only step in pdd, we have approximately pdd — (no dou- 
bly occupied level) for = 1 and poo — (deep level) for 
ed = —5, indicating that Eqs. H17|l and H20fl are valid in 
this bias range. Therefore, we can utilize the deflnition of 
the spin relaxation rate in these equations to account for 
the importance of temperature and bias on the spin-flip 
scattering, as depicted in the insets of Figs. 4(b) and (e). 



III. COUPLED QUANTUM DOTS 

Now we turn to resonant tunneling through a CQD 
with weak coupling between the QDs and the leads, as 
shown in Fig. 1(b). The presumption that the interdot 
hopping is also weak keeps each level of the dots iso- 
lated. Then the superposition of the two levels in differ- 
ent QDs plays a crucial role in tunneling. In order to sim- 
plify our derivation, we consider here the inflnite intradot 
Coulomb repulsion U' and a flnite interdot Coulomb in- 
teraction U, which excludes the state of two electrons in 
the same QD but two electrons can occupy different QDs. 



A. Model Hamiltonian and slave-particle 
representation 

The tunneling Hamiltonian for the CQD is 

(T,cj' k,a 

+ E(^«-4fe.C2.+H.c.), (24) 

where <^n2)a^ Ci{2)cr are creation and annihilation opera- 
tors for a spin-cr electron in the first (second) QD, respec- 
tively, ej [j = 1,2) is the bare-level energy of electrons 
in the jth QD, en2) = i i5j in which 5 is the bare mis- 
match between the two bare levels. The first term in the 
second line denotes the hopping t between the two QDs. 
The other notations are the same with those in Sec. II. 

In the situation interested here, the bare mismatch 5 
should be very small. Otherwise, the quantum coherence 
(the superposition of the two states) has quite tiny effect 
on the tunneling processes. In experiments, this small 
mismatch could be controlled by external time-dependent 
voltages. The available states and the corresponding en- 
ergies for the isolated CQD are: (1) the whole system is 
empty, |0)i|0)2, and the energy is zero; (2) the first QD 
is singly occupied, |cr)i|0)2, and the energy is ei; (3) the 
second QD is singly occupied, |0)i|cr)2, and the energy 
is £2; (4) both of the QDs are singly occupied, |cr)i|(T')2, 
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and the energy is 2ed + U . With the same theoretical 
point of view as in the single QD mentioned in the above 
section, we can decompose the real electron operator Cja 
in these Fock states as 

Cl. = |0)l|0)22(0|l(cT|+^|0)i|a')22('T'|l(a|, (25) 

a' 

C2a = |0)l|0)22(0|l(a|+^|a')l|0)22(fT|l(a'|, (26) 

cr' 

in association with the completeness relation 

|0)l|0)22(0|l(0|+^|(7)i|a')22(CT'|l(^T| 

a, a' 

+ 5I(l^)l|0>22(0|(a| + |0)i|(7)22(^|l(0|) - 1. (27) 

Again, we assign these Dirac brackets as operators: the 
slave-boson operators ~ |0)i|0)2, d^^, = |cr}i|cr')2 
and the pseudo-fermion operators fl^ = |(t)i|0)2, f\„ = 
|0)i|o-)2- Then, Eqs. |(53-(|2IIl can be replaced as 



(28) 
(29) 



+i E(/i-/2. + /L/i J + (2fd + t/) E 

k.a a' 

+ E[^«-4fc-(^^/2. + E fl'd^-^) + H-c.]. (32) 



k.a 



B. The quantum rate equations for the CQD 



Define the retarded (advanced) and lesser (greater) 
Green's functions (GFs) for the CQD G'[j^^'<^>^(i, t') = 
{{Ci^{t)\c]^{t'))Y'^'''^'0^ as usual. Considering Eqs. ^ 
and (|26(l . these GFs can be expressed in terms of the 
slave particles: G^j„ = Geija +E^v" Gdjjaa'a" [* = 2(1) 
if t = 1(2)] with G«,, ^ {{eHt)f.At)\fUt')e{n)) 

and Gdiw." = ((/l.(i)rf.vWl4".(^')/i."(i'))), 

G<i22..'." = {{fUt)d,.,{t)\dl,„{t')f2."it'))); In 
the following derivation, we will use the nondiagonal 
doubly-occupied-related GFs, for example, Gd2ia-a-'a-" = 
{{fL4t)d,'."mlAt')fiAm and G:,2i..v" = 
{{fL{t)d,.'ml„,,{t')f,^„{t'))). Under the weak cou- 
pling assumption and small bare detuning (5, the decou- 
pled GFs of the CQD can be defined in terms of spectrum 
expressions, in the Fourier space, as 



e^e + E(/il/i. + /2./2J + E 4.'^..' = 1- (30) 



And obviously the explicit (anti)communicators of these 
auxiliary particles are: 

ppt — 1 f] , , — () , A , f — /)■ S , 

ee — 1, "criCT2"(jJcr^ — ''(Tier; "cr2 0'2 ; JuyJja' ^"«j"cro- i 



G 



<o 



'^dllaa'a" 
^d22aa'a" 



G 



>o 



= 5a'a"2'Kipddc,'a5[uJ - (cd + U)], 
= Sa'a"2Tripddaa'S[uJ - (cd + U)], 

= -Sa'a"2TTipiia'S[uj - (ed + U)], 



(33) 



edl, = efl = f^^eJ ^ f^J^,^„ ^ d„^,e^ = d,,,„fl = 0. 

(31) 

The density matrix elements are expressed as 
Poo = |0)i|0)22(0|i(0| = ete, pn, = k)i|0)22(0|i((7| = 

/L/l., P22a = |0)l|a)22(cT|l(0| - flf,,, 

Pddaa' = |cr}i|o-')22(cr'|i(CT| = dl^,d^^,, and 

/5i2<7 = 1 0)1 1 (7)22(0 1 1(0- 1 = fLfia- In terms of these slave 
particles operators, the Hamiltonian for the CQD can 
be rewritten as 

^ E ^Vkaclk^C^ka + E •^l^'^^l^ + ^2 E fLf2a 

. I 



Gg<°v"(^) = 0, G>°,(a-) = 0. 

G^j°^(w) = 2mp,.j„6{uj-ed), 
G%l.A^) = -5aa"2^^p■J^J[u-{ed + U)].{M) 

In order to get the quantum rate equations, we use ex- 
actly the same procedure as in the previous section, eval- 
uating the statistical expectations of the rate of time 
change of the density matrix elements pij . After tedious 
but straightforward calculations, eventually we obtain in 
the wide band limit 



POO = dujY,{rLafL{^)G>M + rU^-fL{uj)]G<,,{iu) + rnafR{^)G>^^^^^ 

(35a) 
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Pile 



PT2g 



- ^i^<.'[l-/flM]G^ll.'.."H} + *^(Pl2.-P2l.), 

cr%cr" 

- XI rL^'[l-/i(^)]<^d22<T'aa"('^)}+«^(P21<T-pi2^), 



(35b) 



(35c) 



Pl2o 



^(£2 - ei)pi2. + ^ fdcj{ J2{r„af,MG>,M + r„.[i - UH]G<M} 

V 

a' ,(7" 

- E {rfl.'/fl(^)G',^l..V"(^) + ^rA^ - fR{^)]G'^21aa'AA}} + *t(Pll. - P22.), 



(35d) 



Pdda 



27r 



/ dc.X{ri,/i(a>)G>2,,„,(c.)+rL.[i-/L(^)]G<2,,,,+rfi,/fi(t^)G>i,^^^ 



(35e) 



PddcTcr — 

ZTT 



1 rfc.X{ri^/i(a;)G>2,^,,(c^)+ri.[l-/i(c^)]G<2.^^,+rfl^/H(a;)G>i^^^^ 



Substituting these correlation GFs with their decoupled 
formulations Eqs. and the quantum rate equa- 
tions can be obtained 

(36a) 

Pllo- = rj^^poo + YT Pddaa' - ^LaPllcr 
a' 

~2t3pi2., (36b) 

cr' 

a' 

~ X! ^i<T'P22a + 2t3pi2a, 



Pl2<T = «(e2 - ei)Pl2(T + «<(Pll(T - P22ct) 



(36c) 
(36d) 



Pddaa = r+^pii„ + r+^p22<T - (F^^ + r^^)p<jd„o-, (36e) 

Pddaa = ^RgPlla + ^taP22S ~ (f + ^Ra)PddaS, (36f) 
and along with P00+I]<^(pilcr+P22o-) + I]cr,,T' Pddaa' = 1, 

inwhcihr^^ = V,^af^[f-d) and T^^ = T^„f^{ed+U) have 



r 



(35f) 



the similar prescriptions as in the SQD. In addition, the 
classical parts of the diagonal elements equations have 
the similar interpretations. For example, Eq. (|36bl) for 
the rate of change of the number of the spin-tr electrons 
in the first QD pno- is contributed, noting the fact that 
the first (second) QD do not directly connect to the right 
(left) lead, from four single-particle tunneling processes: 

1) tunneling into the QD with spin-cr electrons from 
the left lead if the QD is initially in the empty state poo; 

2) tunneling out from the QD with spin-cr' electrons F]^^, 
into the right lead if the QD is initially in the doubly oc- 
cupied state Pddaa'] 3) tunneling into the QD with spin-cr' 
electrons F^g.' from the right lead; and 4) tunneling out 
from the QD with spin-cr electrons TJ^^ into the left lead, 
when the QD is initially just in this state pna- Tunnel- 
ing events 1) and 2) increase pna, but events 3) and 4) 
decrease this probability. The final term in Eq. Ij36bl) is 
responsible for coherent effects. The equation Ij36d|l for 
the nondiagonal matrix element pi2a indicates that the 
role of the leads is to provide damping of the quantum 
superposition^ It is also worth noting that the present 
proposed quantum rate equations are reliable for a wide 
range of temperature and external bias voltage, where 
the three major approximations we use are valid. 
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The electric current II flowing from the lead L to the 
QD can be calculated as: 

+ E G,>22..v"(^)] +r,.[l - /,(a;)][G<i,(^) 
+ E G,<22..v"]}- (37) 

ct'.ct" 

Under the weak coupling approximation, it becomes 

Il/R = ^\^l„{pddaa + Pddas) + ^laPlla 
a 

-fl^{p22a + P22a) - ^taPOo]- (38a) 

Similarly, for the current flowing from the lead R we have 

-^-R/e = E {Pddaa + PddSa ) + r^^p22cr 

(J 

+ Piia) - ^RaPaa]- (38b) 

It is easy to prove that, in stationary condition, the cur- 
rent conservation is fulfilled II — —Ir- 

C. Discussions 

In order to simplify the analysis, we only consider spin 
independent tunneling processes in the following discus- 
sions. Two special cases, no doubly occupied state and 
no empty state, are studied. First we assume the inter- 
dot Coulomb interaction U is infinite, whereas only one 
electron can be found inside the system, so Pddaa' — 
and r+, ~ 0. The quantum rate equations (|36b|l - l)36f|l 
simplify to 



It is interesting to compare our results in this situation 
with those of Gurvitz and Prager- 1 for the case of large 
bias voltage between the two leads. For example, the 
large bias voltage determines F^ — 0, F^ = and FJ = 
Tl, T~o = Tr at eV > T. Therefore, the dc current 



becomes 



<2(2 + Tr/2Tl) + (Fi^)74 + (£2 - ei)- 



r, (41) 



which coincides with the result obtained by Gurvitz and 
Pragerii^ It is quite obvious that the finite temperature 
plays a crucial role in the "coherence" tunneling. The 
previous formulations for large bias voltages, however, 
can not provide any information about the temperature 
effects. This is the central improvement of the present 
approach for the coupled quantum systems. 

Secondly, wc consider the deep level situation where 
the bare levels ei and €2 are far below the Fermi level 
but ed + U is just above the Fermi level in equilibrium. 
In this case the CQD is always occupied and poo = 0, 
F~^ ~ 0. Therefore, we have 



pn = 2F^p<irf-2F+pii -2iSpi2, (42a) 
P22 = 2r-p,M-2r+p22 + 2iSpi2, (42b) 
Pdd = r+pu+f+p22-{fl+f^)pdd, (42c) 



P12 = i(e2 -ei)pi2 + ii(pii -P22) - (rj + ^7^)^12, (42d) 

with 2pii -I- 2p22 + 4pdd = 1. After solving the set of 
equations in the steady state, we obtain 



Pll 


= rjpoo - r^pii 


- 2tSpi2, 


(39a) 


Pll 


= K^R + 


i2(F++F+)(F^+F^)/A]/A, (43a) 


P22 


= r^poo - r^p22 


+ 2t5pi2, 


(39b) 


P22 


= [^l^R + 


t2(F++F+)(F^+F^)/A]/A, (43b) 










P12 


- tin^R - 


-F-F+)[ei-e2 + *(r++F+)]/AA, 






\i^l+^R)Pl2, 








(43c) 




)pi2+ii(pii-p22)- 


(39c) 


A 




-2rzr+ + 4F+F+ 














+ F^,+F++F+)(F++F+)/A, (43d) 



with Poo + 2pii + 2p22 = 1. The steady solutions are 

Pll = [rJr^ + <'(r++F+)(F-+F^)/A]/A, (40a) 
P22 = [F-F++t2(F++F+)(F-+F^)/A]/A, (40b) 

P12 =i(F+F^-F^F+)[ei-e2+*i(rZ+r^J]/AA, (40c) 



^ = r^F- + 2F-F+ + 2F+F- 

+t\Tl + F^ + 4F+ + 4F+)(F^ + F^)/A, (40d) 

in which A = (62 - eif + (F^ + T^f/A. The steady 
current is given by Il/e = 2(F^pii - Fjpoo). 



in which A = (e2 — ei)^ -f (f J -I- F^)'^. The dc current is 

lLle = A[flpM-T+p22). 

It is also interesting to consider the situation of large 
bias voltage in the strong interdot Coulomb repulsion C/, 
whereas the Fermi level of the right lead plr lies far below 
Cd + U, but far above the resonance level to satisfy the 
requirement of deep level, meanwhile the Fermi level of 
the left lead p^ is far above -I- [/, so that F^ = 0, 
r+ = and r+ = Fl, F^ = Tr. Finally we obtain 



h/e = - 



2^2(1 -fFi/Ffl) + 2(62-61)2 + (Fi)2 



(44) 
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D. Numerical results 

In this subsection, we perform numerical calculations 
for the tunneling transport through the CQD, by using 
the quantum rate equations (|36|) . in the stationary con- 
dition. We symmetrically add the bias voltage again be- 
tween the source and drain /i^ = —hr = eV/2. 

First we consider the spin independent transport. 
Fig. 5 demonstrates the nonequilibrium occupation num- 
bers in the first and the second QDs, calculated from the 
obtained expectation values of density-matrix elements 

nia = PUa + J2a' Pddaa' and n2a = P22a + J2cr' Pdda'a, 

and the corresponding current versus the discrete level 
for the hopping t — 1.0 between the two QDs at a small 
bias V = 1.0 and a large bias V = 10.0, respectively. We 
find a similar characteristic as in the SQD (Fig. 2): 1) 
The nonzero bias weakens the Coulomb blockade effect; 
2) riio- and n2a- have fractional steps and 3) ni^ ^ 
in nonequilibrium; 4) The conductance has two peaks 
at the resonant points, while the current has three steps 
in the strong nonequilibrium regime. Here we observe 
a higher peak magnitude and a higher step value at the 
deep level regime than those at the no doubly occupied 
level regime, and the maximum step value located at the 
middle "window" of the bare level. More interestingly, 
a opposite behavior has been found when the hopping 
t between two QDs decreases, as shown in Fig. 6(c), in 
which we plot the corresponding results for a small hop- 
ping t — 0.5. Generally, one may expect that increas- 
ing the hopping t can reduce the difference between two 
QDs, and very strong hopping can finally give rise to the 
formation of covalence. In other words, the difference 
between rtio- and 7120- should raising with decreasing the 
hopping t. This is the case as shown in Figs. 6(a) and 
(b), where the occupation numbers are displayed for the 
smaller hopping t — 0.5 in comparison with the results 
of the hopping t = 1.0 in Figs. 5(a) and (b). One can 
note that the occupation number in the second QD even 
experiences a descendance in the middle "window" of the 
bare level for the case of t = 0.5, which expresses an op- 
posite behavior in the case of t = 1.0. This is the reason 
why current-voltage characteristics are different in the 
two cases. 

The effect of the hopping on the tunneling is more 
clearly illustrated in Fig. 7, where we plot the occupation 
numbers, the current, and the differential conductance as 
a function of the bias for different hoppings t in the no 
doubly occupied level = 1 and the deep level = —5. 
In both cases, we have two peaks in the differential con- 
ductance corresponding to the two steps in the current. 
More importantly, we find that the current declines in 
the second step and consequently the negative differen- 
tial conductance (NDC) appears in the according biases 
when the hopping t < 1.0. We can explain appearance of 
the NDC by variations of the occupation numbers with 
the bias, as shown in Figs. 7(a) and (d) for the hoppings 
t = 0.2 (thick lines) and t = 1.0 (thin lines). Consid- 
ering the fact that we apply the bias symmetrically and 



/i/i = —eV/2 < 0, the current flowing from the right 
lead is dominated for the case = 1 by the process: 
tunneling out from the second QD into the right lead. 
According to Eq. (|38b|l . we have Ir/c « J2a^R'^''^'^'^- 
is obvious from Fig. 7(a) that the rising second step in 
7i2(T for the case oi t = 1.0 (thin dashed curve) indicates 
the rising step in the current, whereas the decline second 
step for the case oi t = 0.2 (thick dashed curve) implies 
the NDC. In the other case — —5, the current flow- 
ing from the left lead is ruled by the tunneling process 
into the first QD from the left lead, being approximately 
/i/e = J2cr rLCT[l — 'T^iT Apparently, the variations 

of nia denoted by the solid hues in Fig. 7(d) provide 
interpretations for the current-voltage characteristic in 
Fig. 7(e) and the NDC in Fig. 7(f). Therefore, it can be 
addressed that open of a new channel provides negative 
contribution to the current in the case of t < 1.0. 

An interesting question is what happens for the tun- 
neling current and the NDC when the interdot Coulomb 
interaction U weakens or strengthens. We show this in 
Fig. 8, where current vs bias is presented for various cor- 
relation parameters from J7 = to 00 in the cases t — 1.0 
(thick lines) and t = 0.2 (thin lines). For U 00, the 
current has only one step with increasing bias in both 
cases of = 1 (a) and = —5, because no new channel 
is available due to the extremely strong charging effect. 
For the finite interdot Coulomb correlation, however, the 
applied bias can overcome the Coulomb blockade effect 
and open a new channel for tunneling at the correspond- 
ing threshold value of voltage. For t < 1.0, this new 
channel induces a peak in the current. This peak be- 
comes narrower with declining value of U, but its height 
remains unchanged if U is not too small. At sufficiently 
small values of U, as shown in the inset of Fig. 8(a), 
height of the peak in current decrease, even vanishes fi- 
nally when U = 0. So we can claim that the interdot 
Coulomb interaction U = and 00 leads to the single 
peak in the differential conductance, but the finite val- 
ues result in double peaks, and even the NDC in the case 
of t < 1.0. 

The temperature effect is also shown in Fig. 8 for 
[/ = 4. Increasing temperature smoothes I-V curve, but 
remains the NDC unchanged. 

Now we study the spin dependent tunneling through 
the CQD connected to two ferromagnetic leads. Figs. 9 
and 10 depict the occupation numbers in the two QDs 
and the current in both P and AP configurations for 
6(1 = 1 and erf — —5, respectively. We find, besides 
analogous behaviors with the spin-independent tunnel- 
ing, that: 1) rill ^ riii even in both alignments; 2) 
n2] — n2i in the P configuration but n2i 7^ ri2| in the 
AP configuration; 3) variations of the current flowing in 
different magnetic configurations are very sensitive to the 
value of the hopping between two QDs, which leads to 
4) the negative TMR for the sufficiently small hopping 
t = 0.2, as exhibited in Figs. 9(c) and 10(c) at certain 
voltages. 
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IV. CONCLUSION 

In this paper, we have systematically derived the quan- 
tum rate equations for sequential tunneling from NGF, 
and then utilized them to investigate quantum coherent 
transport in a single QD with weak spin-flip scattering 
and weakly coupled QDs systems taking the intradot and 
interdot Coulomb interactions into account. In these sys- 
tems, the superposition between different states plays a 
vitally important role in coherent tunneling processes. 
Directly, a kind of quantum oscillations in mesoscopic 
systems is due to this superposition effect. Now, it is 
believed that the master equations or the modified quan- 
tum rate equations, which are actually equations of mo- 
tion of density submatrix for diagonal and nondiagonal 
elements, provide a successful tool to study this phe- 
nomenon, and even allow an analytical description. 

For this purpose, we have generalized the slave particle 
technique, which is developed previously in the single- 
site space and successfully applied to study the strongly 
correlated systems, into the two-site space. Based on 
this theoretical approach and the correct quantization of 
these artificially introduced operators, previously well- 
developed NGF for noninteracting systems has been used 
to construct the quantum rate equations when only three 
assumptions are made: first, the coupling between the 
central region and the leads must be weak; second, the 
couplings between the subsystems are also weak, for ex- 
ample, weak spin-flip scattering in SQD and weak inter- 
dot hopping in CQD; third, the wide band limit. The first 
condition makes it valid that we can keep only the lowest 
order terms in \V\'^ in the expansions of the equations of 
motion. It also renders the central region approximately 
a qusi-equilibrium isolated system, which facilitates the 
"localized" energy spectrum expressions for the corre- 
lation GFs of every subsystem in the central region in 
combination with the second presumption. These ap- 
proximations not withstanding, our approach is appro- 
priate for a wide range of temperature and external bias 
voltage, and incorporation of the charging effect. Finally, 
it should be pointed out that our derivation is equivalent 
to the lowest-order gradient expansion techniqucn^ 

Employing this approach, we have studied in detail the 
coherent tunneling through a SQD and a CQD systems. 
We have given some analytic expressions for steady- 
state transport in two special cases: doubly-occupied 
prohibited state and deep level in large intra- or inter- 
dot Coulomb repulsion. Furthermore, we have compared 
some of our results with previously obtained results in the 
literature. For example, for resonant tunneling through 
a SQD with spin-flip scattering, our approach provides a 
quantum correction to the classical results. When there is 
no spin-flip scattering, our rate equations reduce exactly 
to the classical results as Glazman and Matveev)^ and 
Beenakkeri^ In the case of resonant tunneling through a 
CQD, our results are in perfect agreement with the pre- 
vious analysis proposed by Gurvitz and Prager^5. under 
the limitation of zero temperature and large bias voltage. 



In addition, we have performed numerical simulations 
for variations of occupation numbers and the current with 
increasing bias voltage and varying the discrete level in 
QD. We summarize the main common features as follows: 
1) Occupation numbers have frational steps in nonequi- 
librium, implying that the Coulomb blockade effect is 
partially overcome by applying bias voltage; and corre- 
spondingly 2) the current-voltage characteristic displays 
two steps, giving rise to double peaks in the differen- 
tial conductance. Especially, our calculations manifest 
the importance of temperature and bias on the spin-flip 
transitions in the SQD. For the CQD, a possible NDC 
can be reached if the interdot Coulomb interaction is fi- 
nite and the hopping between two QDs is small t < 1.0. 
Besides, the TMR becomes negative in nonequilibrium 
for the CQD connected to two ferromagnetic leads if the 
hopping t is sufficiently small. 
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APPENDIX: DERIVATION OF QUANTUM RATE 
EQUATIONS FOR SQD 

In the Appendix, we present a detail derivation of 
Eqs. Hll() . In the following, we take Eq. (jllb|l as an 
example. The statistical expectation of Eq. (|10b|l gives 

+iR%PaS ~ iRaPsa- (A.l) 

According to Langreth's operational rules'^^, those hybrid 
correlation GFs are given by 

+G< ^(t, h)v;,,g'^,,,, {h,t% (A.2a) 

+G<^At^h)V;,,g'^^„,{t,,t% (A.2b) 

G:ik.'.eaM = '5^^' J dt,[g:^,,,{t,t,)V,^,Gf,,{h,t') 
+9^k.'it^t')Vv^'G:,Ati,t% (A.2c) 

(A.2d) 
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Substituting Eqs. HA.2|) into Eq. HA.1|I and taking the 
Fourier transformation, pcro- can be expressed as 

r/.k 

+G<,M(g^fe,M-g;fe,M)] 

+iRgPaa - iRcrPaa, (A. 3) 

where grfkai^) are the Fourier transform of the exact GFs 
in the 77th lead without the couphng to the central region. 



In the wide band limit, one has 

J2\V,af9^,A^) - ^T,,,Mu;), (A.4a) 

k 

ElVlV/c^H - -*r,,[i-/„H]. (A.4b) 

k 

Substituting the GFs (|A.4p into Eq. (|A.3|I and employing 
C -G" = G> -G<, Eq. fTTE)l can be reached. Analo- 
gously, we can derive other equations in Hll|) . 
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FIG. 1: Schematic diagrams for the resonant tunnehng 
through (a) a single interacting QD and (b) a coherently cou- 
pled QDs. 
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FIG. 2: Noncquilibrium occupation numbers n|, n|, and pdd 
(a,b), and tunneling current (c) vs the bare level of the SQD 
with no spin-flip scattering for both magnetization configu- 
rations, (a) is plotted at a small bias V = 1.0 and (b) is at 
a large bias V = 10.0. The thick lines are plotted for the 
AP configuration, and the thin curves are for the P configu- 
ration. The equilibrium occupation rmrnbcrs are depicted by 
the thin lines in (b). Other parameters are: U = 4, T = 0.1, 
and p = 0.5. 
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FIG. 3: Occupation numbers n|, rij (a,d), pdd (b,e), and 
current (c,f) vs the bias voltage, calculated for no spin-flip 
processes and different temperatures T = 0.1, 0.5, and 1.0. 
(a)-(c) are plotted for ea = 1, (d)-(f) for = —5. The insets 
in (c) and (f): the corresponding TMR vs the bias voltage. 
Other parameters are as in Fig. 2. 
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FIG. 4: Occupation rmmbors n^, ni (a,d), pdd (b,e), and 
current (c,f) vs the bias voltage calculated for the AP con- 
figuration with the spin-flip transition R = 1 and different 
temperatures T = 0.1, 0.5, and 1.0. (a)-(c) are plotted for 
€d = 1, (d)-(f) for td = —5. For comparison, the respective 
results without the spin-flip transition are also plotted as thin 
lines. The insets in (b) and (c): the temperature and bias de- 
pendence of the spin relaxation rate. Other parameters are 
as in Fig. 2. 
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FIG. 5: Noncquilibrium occupation numbers nia-, n2cr, pdd 
(a,b), and current (c) vs the bare level of the CQD. (a) is 
plotted at a small bias V — 1.0 and (b) is at a large bias V = 
10.0. The equilibrium occupation numbers are also depicted 
by the thin lines in (b) and (e). Other parameters are: J7 = 4, 
T = 0.1, and t = 1.0. 
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FIG. 6: The same as Fig. 5 except for t = 0.5. 
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FIG. 7: Occupation numbers nia-, n2a, Pdd (a,d), current 
(b,e), and the differential conductance (c,f) vs the bias volt- 
age, calculated for different hopping t. (a)-(c) are plotted for 
ed = 1, (d)-(f) for ed = —5. Other parameters are as in Fig. 5. 
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FIG. 8: The current-voltage characteristics, calculated for in- 
dicated values of the interdot Coulomb interaction U and dif- 
ferent hoppings t = 0.2 (thin lines) and 1.0 (thick lines), (a) 
is plotted for = 1, (b) is for = —5. Inset in (a): I-V 
curves for small U . 
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FIG. 9: Occupation numbers riij, ni|, n-z], and n2i in the P 
configuration (a) and the AP configuration (b) for t — 0.2, 
current (c), and TMR (d) for t = 0.2 and 0.5 vs the bias 
voltage. Other parameters are = 1, T = 0.1, and p = 0.5. 



1.0 
0.8 

c 

•2 0.6 

3 0.4 

o 

o 

O 0.2 



f — 



0.0 - 



0.2 



0.1 - 



0.0 



(a) . 



AP 
' 

V." 



r — ^ 



(c) 



AP 



1=0.2 

0.5 



AP 

_l 



(b) 



_L 



T=0.1 



10 
V 



15 



10 
V 



0.3 

(diJ 

- 0.2 



- 0.1 
~- 0.0 
J -0.1 
-0.2 



DC 



15 20 



FIG. 10: The same as Fig. 9 but for the case of td 



